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Abstract We use recently introduced notion of stochastic soliton lattice for quantitative 
description of soliton turbulence. We consider the stochastic soliton lattice on a special 
band-gap scaling of the spectral surface of genus so that the integrated density of states 
remains finite as A^ — ^> oo (thermodynamic type limit). We prove existence of the limiting 
stationary ergodic process and associate it with the soliton turbulence. The phase space of 
the soliton turbulence is a one- dimensional space with the random Poisson measure. The 
zero density limit of the soliton turbulence coincides with the Frish - Lloyd potential of the 
quantum theory of disordered systems. 

1 Introduction 

The idea of introducing the statistical description into the soliton theory has its origin in 
the work of Zakharov on kinetic equation for the rarefied soliton gas. In this work, a 
new object, an infinite sequence of KdV solitons on the x-axis was considered. An accurate 
limit for the A^- soliton solution as A^ ^ cxd is , however, very nontrivial and has been 
investigated so far only for the case of a very special phase distribution ||^ which yields 
the reflectionless potential decaying at one of infinities x ±00. This kind of infinite 
- soliton solutions does not provide directly the spatial uniformity which is necessary for 
description of the soliton gas in thermodynamically equilibrium state. 

Another view on the problem was proposed by Lax ||^ who noticed that 'the weak limits 
of oscillatory sequences of dispersive compressible flows show a remarkable resemblance to 
ensemble averages of classical turbulence theory'. To constitute such a deterministic analogue 
of turbulence one has to construct the correspondent ensemble of flows and provide it with 
the physically plausible measure which form a stochastic process. It is clear from the very 



beginning that for description of the 'sohton turbulence' we are interested in the stationary 
ergodic processes. 

The consistent way of apphcation of stochastic description to the sohton theory was pro- 
posed in p, where the notion of stochastic sohton lattice (SSL) was introduced. The 
basic idea of such a description lies in the fact that the finite-gap solutions of completely inte- 
grable equations are almost periodic functions and posess, therefore, their natural stochastic 
structure determined by the compact shift group with the uniform (Haar) measure. 

In this work, we investigate A^-gap stochastic soliton lattices q{x) = vn^x) of the KdV 
equation on a special band-gap scaling of the spectral Riemann surface of the complex 
parameter E when the width of gaps is ~ 0{1/N) while the bands are exponentially narrow 
~ exp(— A^), ^ 1. Along with this, the small parameter 1/A^ in our consideration does 
not appear explicitely in the Schrodinger equation {—dl^ + q{x))(f) = Ecf). 

The basic results of the present work are: 

• the process C,{x) = lim unIx) exists. It is an ergodic stationary random process 

each realization of which satisfies the KdV equation. The phase space of ({x) is a 
one-dimensional space with random Poisson measure. 

• the considered limit is of a thermodynamic type (the integrated density of states re- 
mains finite as A^ — s> oo for the chosen scaling of the band-gap stucture). 



band 
gap 

interpreted soliton turbulence. 



as lim the thermodynamic limit of stochastic soliton lattices C(x) can be 



• we calculate the rotation number and some ensemble averages for the soliton turbu- 
lence. 

• we show that the zero-density limit of the soliton turbulence yields the Frish - Lloyd 
potential of the quantum theory of disordered systems p, . 

2 Basic Preliminaries 

The soliton lattice (SL) (A^-gap potential of the Schrodinger equation) is given by the known 
Its - Matveev formula (see for instance [0) 

u^{x;r) = C(r) - 29^ loge[z(x)|i?(r)] , (1) 
x e R, r = (ri, . . . , rajv+i), ri < ra < . . . < rgTv < r2N+i ■ 

Here 

&[z\B{r)] = ^exp{7ri[2(m,z) + (m,fim)]}, (2) 

m 

m = (mi,...,m^) G Z^, z G 
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is the Jacobi theta-function of the hyperelhptic Riemann surface of genus 

2N+1 

/i2= n {E~r,) = R^N+i{E,v) (3) 

with cuts along the bands. 

The Riemann matrix -B(r) and the constant C(r) are expressed in terms of the basis 
holomorphic differentials ipj : 

. N N . 

B,, = fij,, C{v) = Y.r,-2j2fEi:,, (4) 

Here 

N-l j^k 

fc=0 JR2N+l{r,E) 



and dependence ajk{r) is given by the normalization 

= 6jk , (6) 



where the a-cycles surround the bands clockwise, and the f3- cycles are canonically conjugated 
to a's such that the contour Pj starts from the j-th cut, then goes to +oo and returns on 
the lower cut. 

The imaginary phases z(x) are given by the formula 

z = —2isiN~ix + d , (7) 

where d is the initial imaginary phase vector. 
Using the substitution 



27r 

we rewrite ( |l|) in the form with real phases y |T2[ 



- = ^By (8) 



UN{x\r) = UN{yi{x), . . .,yN{x)\r) , (9) 



yj{x) = kjX + fj , fj{mod27i) , (10) 

k = -4mB-^aiN^i . (11) 

Here k = k(r) = {ki, . . . , k^) is the wave number vector, and f = (/i, . . . , /at) is the initial 
(angle) phase vector, — tt < fj < tt. Note that in (||) 

MAr(?/l, . . . , ?/j + 27r, . . . , ?/Ar|r) = Ujv(yi, . . . , . . . , ?/Ar|r) , (12) 
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that is UN{x\r) is A^-quasiperiodic in x |T^, (hereafter we consider only incommensurate 



kj,j = 1,...,N). As for any almost periodic function, we have for UN{x\r) the Fourier 
representation 

UM{x\r) = Y^Cje"^'^^+''^\ (13) 
j 

where Cj , Ij , hj are real, —it < hj < tt. {Ij G M, where M is a frequency-module [§] and hj 
are the integer linear combinations of fj{mod27T)). 

The KdV evolution of (^) is isospectral and is described by the linear motion of the 
phases on the Jacobian: 

MAr(x,t|r) = {yi{x,t), . . . ,?/Ar(x,t)|r) , (14) 
yj{x, t) = kjX + ujjt + fj , (15) 

where the frequency vector 

2Af+l 

uj = u{y) = {ui,. . . ,un) = -^TiiB~'^{a.N-i rj + 2&N-2)- (16) 

i=i 

3 Stochastic Soliton Lattices 

It is well known that any almost periodic function generates the stochastic stationary process 
(see i, i). 

Definition 1.0, The stochastic process generated by SL u^{x\y) we call Stochastic 
Soliton Lattice (SSL) and denote as i/Ar(x|r). 

The general construction of z/7v(x|r) admits a very simple and clear description. The 
realization set Vl of it consists of functions (^, ([T0|) where f G Tor^ ; Tor'^ is A^- dimensional 
torus {—71 , vr]^. The probability measure dfi is the uniform (Lebesque) measure on the torus. 
It corresponds to the description of z/iv(x|r) following from(|j),(|TO|): 

z/Ar(x|r) = un{- ■ ■ kjX + (pj . . . |r) (17) 

where 0i, . . . , 0Ar are independent random values uniformly distributed on (— vr , vr], that is 
= (01, . . . , 0Ar) is uniformly distributed on Tor^ . As kj are incommensurate then vn{x\v) 
is an ergodic process . 



As h'N{'x\Y) is the stationary process then it has the Stone - Kolmogorov spectral decom- 
position [|l5l; due to the ergodicity this decomposition has the form (cf. (|13])) 



MA^) = 11^3^'^''''^''\ (18) 



where 6j are uniformly distributed on (— vr , tt] noncorrelated random values |T7 
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The well known formula (Bochner - Khintchin) [^] gives us the covariance function K{h) 
of the stationary process z/7v(a;|r): 

K{h) = (z>^(x|r) ■ On{x + h\r)) = ^ |c,f e^'^^ . (19) 

j 

Here ^(x) = ^(x) — (^) is the centered process. 

Theorem. Consider the KdV equation ut — 6uUx + Uxxx = as the equation describing 
the evolution in the phase space of stationary processes. Let the initial data has the form of 
the SSL: 



u{x,0) = UM{x\r) . (20) 

Then the solution of the KdV is 

u{x,t) = ujy^xlr) = u{x,0) , (21) 

that is u{x,0) = z/Ar(x|r) is a stationary point. 

Proof . The evolution of realizations (^ ,(p!OD is described by ([T^) ,([T5|). For any moment 
t one can introduce the new 'initial phase' f* = Ujt + fj which is also uniformly distributed 
on Tor^ . Therefore, the KdV evolution changes neither realization set Vt nor the probability 
measure. Q.E.D. 

Now we present some expressions for the ensemble averages which will be needed in the 
future. If Q{f) is an arbitrary smooth function then the mean value Q(0 of stochastic 
process ^ is 

{Q{0) = JdMf), fen. (22) 

n 

(We recall that ^(x) = {^7 = {/(x)} ,B,fi}, and B is some cr-algebra of measurable Borel 
sets of n). 



For z/7v(x|r) we have 12 



_^ IT TT ^ L 

((5(z/Ar(x|r)) = j^—y^ I "' I ^'^1 • • • #AfQ(i^7v(a;|r)) = lim^ - J (5(z/Ar(x|r)) , (23) 



— 7r —IT 



which constitutes ergodicity of SSL. Direct calculation using ([|),(0) gives surprisingly simple 
formulas for two first moments of the value z/7v(x|r) Namely, 

{uN{x\r)) = C{r) - -{k, Bk) , (z/^(x|r)) = -^^{uj, Bk) , (24) 

and the vectors k(r),co'(r) are given by ( pH]) , ([T6|). It should be noted that expressions 
are obtained for the particular choice of the canonical basis of cycles (see (^)). However, 
namely this normalizatiton is preferrable for our consideration. 
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4 Rotation Number in Stochastic Soliton Lattices 

Consider the Schrodinger equation with an almost periodic potential q{x) 

{-dl, + q{x))<l) = E<f), xeK. (25) 

The potential q{x) has an important characteristics, the rotation number, which is defined 
for real E as [JM] 

a(E) = lim - arg(0'(x, E) + i^(x, E)) . (26) 

We will use also the integrated density of states M {E) which is connected with the rotation 
number by a simple relation (here we improve some inaccuracies in j^, P): 



U{E) = -a{E) . 



TC 



If q{x) is N-gap potential , q{x) = UN{x\r), then 



(27) 



ttNiE) = Re J dp{E') 



[28) 



where dp{E) is the quasimomentum differential [Q, |]TT 
E^ + foiv-i^"^-' + ...+_bo 

^3 



dp{E) 



R2N+i{E; r) 



fdp{E) = 0, j = l,...,iV. (29) 



The integrals of dp over the a - cycles are known [|T2[, |TT| to give the components of the 
wave number vector ([TTD 



^ dp{E) = kj . 



(30) 



Then, 



1 M{E) 



if E G gapA^ 



(31) 



1 M{E) E 

- ^3+ / ME') ii E e bandM 

2 j=i J 



where M{E) is the number of the band nearest to E from the left M{E) < N , M(0) = 
N, M(-l) = 0. 

We note that the formula ( |3TD gives an effectivization for the case of finite-gap potentials 
of the statement in |^ that the values of the function 2a{E) ii E E {gap} belong to the 
frequency module of the almost periodic potential. 
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As the SL (I), ([T0|) is the quasiperiodic function both in x and t one can also formally 
introduce the temporal rotation number for the SSL by 

E 

a'j^{E) =Re J dq{E') (32) 

— oo 

where dq{E) is the quasienergy differential fT^'O {dq{E) = (E^+^ + cnE^ + ... + 

Co)dE/ J R2N+i{E] r), / dq{E) = , / dq{E) = Uj, where Uj are the frequences (pJ])). 

ft 

Since the value of a^iE) is the same for any realization from the set Q from Def.l then 
it is the spectral characteristics of the whole SSL z/jv(x|r) (deterministic characteristics of 
the stochastic process). Due to the Theorem 1 the rotation number does not change under 
the KdV evolution. 

We introduce also the full integrated density of states in the SSL A which is a number: 

1 ^ 

X^Ar^{oo) = —J2k,. (33) 
i=l 

One can see that the full integrated density of states has in the SSL the natural meaning of 
the mean number of waves per unit length. 



5 Thermodynamic Limit of Stochastic SoHton Lattices 

Let the nontrivial (band-gap) part of the spectrum lies in the interval (—1,0) of the real 
£'-axis: aAr(— 1) = 0, a^iO) = ttX. For iV 1 we consider the following scaling of the 
band-gap structure 

gaps(E) ~ ^ bands(E) ~ exp (-iV) > 1 , E e (-1, 0) . (34) 

Using (|TT|) one can show that the scaling ( |5^ ) implies the following behavior for the wave 
numbers: kj ~ which provides finiteness of the rotation number a^iE) (and of the 

integrated density of states Nn{E)) as N ^ oo: 

I M{E) 

a{E) = lim a7v(^) = lim - V A;. < cx) , M <N. (35) 

N^oo N^oo 2 — ^ 



Due to this property we call the scaling ([3^ ) and the correspondent limit as — oo the 
thermodynamic ones. 

It should be noted that the thermodynamic scaling (0) appears when one considers the 
quasiclassical asymptotics of the spectrum for periodic potentials [0, [|^, [^. In contrast 



to indicated works, however, the small parameter 1/A^ in our consideration is contained only 
in the band-gap structure of the spectrum and does not appear explicitely in the Schrodinger 
equation (pSj). As a result, we will show that the limit for the finite-gap potentials on the 
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thermodynamic scaling exists as N ^ oo in a strong sense. We emphasize the difference 
of this point from the known Lax - Levermore - Venakides (LLV) approach (see 
pop to the hmiting passage hmTv^oo^Ar where due to rescahng x — > Nx the above limit 
exists only in a weak sense. The LLV consideration provides information about the slow 
modulations {macro structure) appearing in the A^-soliton (or A^-gap) KdV solution under 
the time evolution. We, on the contrary, are interested in the limiting micro structure of 
nonmodulated A^-phase wave as iV — > oo. Futhermore, we are going to show that the 
thermodynamic limit exists for the whole stochastic process (SSL) which implies not only 
compactness of the realization set but also convergence of the probability measure. 

Now we describe the thermodynamic scaling more in detail. Following we introduce 
the lattice of points 

1 ^ r/i > r/2 > • • • > ^jv ~ , (36) 
where ^ 

-Tfj = 2 ir2j~i + r2j) 

are centers of bands. 

We define two continuous functions on the lattice: 

1. The normalized density of bands: ip{ri)drj k, ^r o — an s m (E,E+dE) 

1 1 } 

^iVj - Vj+i) N 7 

2. The normalized logarithmic band width 

liVj) = log(r2, - r2,_i) + 0(1) . (38) 

Now one can easily represent the relationship (^) for the thermodynamic limit of the 
rotation number in a continuum form: 

2da{-rj^) = <^{r])k{ri)dri . (39) 

where k{ri) is the continu of the wave number vector ([ll|). 

For the temporal rotation number a\f{E) ( p2D in the thermodynamic limit we have the 
equation which is analogous to (pOf): 



2da\-r{^) = ip{ri)uj{ri)dri . (40) 

where w{ri) is the continuum limit for the frequency vector (p!6|) . 

To find a continuum limit for the wave number and the frequency vectors given by (jTT|) 
and ( p^ one has to know the continuum limits of the period matrix B (^) and of the vectors 
a^v-i and ajv_2 from (^. It is clear that all quantities given on the Riemann surface can be 
computed now in terms of two functions (densities) ip{r]) and 7(77) completely describing the 
band - gap structure. 
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2N+1 



In particular, the period matrix (H) and the vectors aAr_i and (aAr_i + 2aAr_2) 



take the form ||20ll: 



5^fi(r/,/x) = — lo^ 



TT 



where 5{x) is the Dirac delta-function 



JL 

'2n 



{log 


T)- IJ, 








2N+1 




1 E 







(41) 



TT 



(42) 



Now we represent the general relationship (|TTD between the wave number vector k and 
the vector aAr_i in the form suitable for the limiting transition 



hB = —AniaN^i 



(43) 



Then applying (pOf), (^T]), (^) to (|i3| ) we arrive at the integral equation for the rotation 
number in the thermodynamic limit : 



vr 



1 

/log 


r]-fi 





r] + fi 



da^—jj^) + 'j{r])da{—r]'^) = —r]Lp(r])dr] 



(44) 



Similarly, for the thermodynamic limit of the temporal rotation number one gets with 
the aid of (0), (|1|), (0) and (|1§: 



1 

— ( 

TT ' 



1 

/log 


?7-/x 


da* 










da\-i/) + '^{'q)da\-7f) = Aifif{7])dri . 



(45) 



Thus, we have established the existence of the thermodynamic limit for the rotation 
number in SSL (as a matter of fact, the rotation number, as well as in the finite-gap case, is 
a deterministic function). 

The rotation number is known to be an imaginary part of the Floquet exponent w{E) 
(for the finite-gap potentials w{E) = i J dp{E)) P). Then we have with the aid of Herglotz 
formula: ^ 

a{z)dz 



w{E) 



E-z 



ImE > . 



But following Kotani P2[, if the Floquet exponent satisfies the conditions 

Rew{E + iO)=0 a.e. on [0, +oo) 



/ -Rew{E + iO)da{E) = , 



(46) 



(47) 
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then there exists an ergodic stationary random process whose Floquet exponent is w{E). One 
can see that for SSL conditions (^) are satisfied due to the properties of the quasimomentum 
(p9|). Therefore the thermodynamic hmit of SSL exists. 

Remark. It can be shown that the existence of the hmit for the rotation number 
guarantees the existence of the hmit for the Weyl function. Then the compactness of the 
reahzation set can be deduced from the results of Marchenko pSf . 

As band/ gap — >■ in the thermodynamic hmit (see (^^), then it is natural to call the 
thermodynamic limit ({x) of the SSL unIx) a soliton turbulence 

C{x) = lim un{x) . (48) 

We believe that the spectrum of a typical realization of the soliton turbulence is pure point 
on (—1, 0) but exact proof of this conjecture is missed yet. 

The ensemble averages (moments) (|22D are the important characteristics of the soliton 
turbulence. To compute three first moments we make use of the expressions (plf) obtained 
for the finite-gap SSL. These expressions admit direct limiting passage as ^ oo in the 
thermodynamic scaling (|34D . 

Using results of it is not hard to show that the constant C(r) (^) in the Its-Matveev 
formula (|l]) vanishes in the thermodynamic limit. Then, substituting (|39D, ( ^OD into (^) we 
have for the two first moments in soliton turbulence: 



(() = -2A J da{E) = -2A , (49) 

-1 

4A r 4A_ 



{C) = YJ Eda{E) = -^E. (50) 
-1 

Here 

, 

da{E) = , lda = l, f{E) = J f{E)da{E) . 



Thus, we have expressed the ensemble averages in soliton turbulence through the averages 
over the spectrum with the measure da{E). One should also note that the Floquet exponent 
w{E) (^61) provides the averaged Kruskal integrals as the coefficients in the decomposition 
as E ^ oo. 



6 Poissonic Properties of Soliton Turbulence 

We study the phase space of the soliton turbulence. We recall that in the finite-gap SSL 
the phase space is the - dimensional torus with the uniform (Lebesque) measure. For 
performing the thermodynamic limit it is convenient to introduce the linear phases instead 
of the angle phases in the SSL (\[7\} 



10 



where Ij [j = 1, . . . , N) are independent random values uniformly distributed on {—f-,f-] 
respectively. 

The further consideration is close to standard procedure of the thermodynamic limiting 
passage in the ergodic theory of ideal gas (see for instance fl^). We now turn to a 



new space with the aid of factorization of the torus by interchange group Sn- The sets of 
points Ij serve now as the points of the obtained space Qn = Tor^ /Sn- The measurable sets 
in the factorized space consist of realizations having 0,1, ... ,N linear phases in the interval 
A of the X axis. 

We introduce the random value = X(o,i)(^i) which is the number of hitting of the 
particular phase Ij into the fixed interval (0,1) C R (we suppose that ir/kj > 1). The 
variable C,j takes two values: 1 and with the probabilities Pj{l) = kj/2'K and Pj{0) = qj = 
1 — Pj = 1 — kj /2tt . The generating function (pj{z) for is 

= - Pj) + zpj (52) 

N 

The sum C,^^' = J2 number of hitting of all linear phases into (0, 1). As C,j are 

independent random values, then the generating function for C,^^^ has the form 

y.(^)(.) = n vA^) = 11(1 + - i)p,) = n(i + ^^-i^) . (53) 

j = l j = l j = l 

On the thermodynamic scaling ( P^ kj = 0{N^^). Then taking the logarithm of (^) we 
obtain 

TV 



1 

ln^W(z) = (^-l)-5:A:, + 0(Ar-^). 
Therefore taking into account (|33D one has 



y.(~)(z) = exp{(^ - 1)A} = E ^" ■ (54) 



n=0 



W. 



Thus, the right-hand part of ( ^4]) is the generating function for the Poisson distribution 
with the parameter A ( which is the full integrated density of states in our case). Therefore, 
the limiting measure in the configurational space is the random Poisson measure (Poissonic 
white noise). 

We note that the Poissonic white noise can be described as the random collection of the 
points {Ij} (the linear phases in our case) on the x - axis such that the distances Sk between 
them are independent random values distributed exponentially with the density 

/(s) = Aexp(-As). (55) 
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We also note that the Poisson parameter A according to (^) can be expressed through the 
first moment in the soliton turbulence: 




The value A can be then interpreted as a density of the soliton turbulence. 



7 The Frish — Lloyd Potential as a Zero-Density Limit 
of the Soliton Turbulence 

We consider the soliton turbulence with the small integrated density of states. We introduce 
a small parameter e <^ 1 and make a renormalization: 



a{E) = ea{E) 



X = eX. 



(56) 



This implies the following transformations for the functions ip{ri) and 7(77) (pTf), (|38| ) 
characterizing the Riemann surface in the thermodynamic scaling 



7(r/)=0(l). 



(57) 



We choose 7(77) = ct] which corresponds to imposing the periodicity condition : \/kj 
2TTc-\eN-^). Then, it follows from (^), (||) that A = 1/c. 
We also make rescaling of the spatial variables 



X 

X = — 
e 



I 



(5^ 



Consider now the SSL with a small integrated density of states on the thermodynamic 
scaling (|5^. We make use of the Its - Matveev formula (|I|): 



where according to (|41|), (|4^ ) 



-2iaN^ij{x — Ij) 



\B] 



(59) 



R. . ~ 



Bij ^ lo^ 



Vi - Vj 



Vi + Vj 



(60) 



(61) 



An analysis of the explicit expression for the theta- function (H) provided (|60|) , (|6lD shows 
that on the set of realizations excluding the set of a small measure one can neglect the 
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contribution of the off-diagonal part of the period matrix B into the solution ( ^9]) (one 
suppose the inequality 



rurij , > 1 . (62) 



j=i ^ i^j=i 



Vi - Vj 



Vi + Vj 



to be satisfied outside of the indicated set of a small measure). 

Then asymptotically in A^, e, {N 3> 1 e -C 1 , Ne ^ 1 ) the following factorization is 
valid : 

e^^ne(r^,^^|r,), (63) 

j=i Tre: 

where Qi^y) is the one-dimensional theta-function (G3 in standard literature (see for ex. 

H)), 

K(l — rrij) 1 , c , , , 

^' = -lfKr=="- ™.«l--exp(--,,)«l. (64) 

Here K{m) is the complete elliptic integral of the first kind, nij are the ellipticity parameters 
(do not confuse with the multiintegers in (H), (|62D ); mj = 1 corresponds to the soliton limit. 

Consider the rescaled SSL of a small density i'N{x)/e'^. Then it follows from (|59D , (|63D , 
(0) that 

(;^\ N sech"^ i Vi^^—^ I 

-2 V^- (65) 



£2 - 5- 



Then passing to a limit 



AT ^ 00 , eN 00 (66) 



one obtains 



^-(5) = lini ^ = -4 E a?5(£ - /,) , (67) 



where, according to the results of Sec. 5, L are random points on the The distances 

Sk between these points are independent random values distributed with the density 

/(5) = iexp(--). (68) 
c c 

The amplitudes aj are independent among themselves and independent of Ij random values 
distributed with the density (p{a). The distribution ( p7D is so-called complex Poisson white 
noise; it is called the Frish - Lloyd potential in the quantum theory of disordered systems 
i, [0. 
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